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Abstract 


In this paper, a genetic-algorithm-based artificial neural network (GAANN) model radioactivity prediction 


is proposed, which is verified by measuring results from Long Range Alpha Detector (LRAD). GAANN can integrate 


capabilities of approximation of Artificial Neural Networks (ANN) and of global optimization of Genetic Algorithms 


(GA) so that the hybrid model can enhance capability of generalization and prediction accuracy, theoretically. With 


this model, both the number of hidden nodes and connection weights matrix in ANN are optimized using genetic 


operation. The real data sets are applied to the introduced method and the results are discussed and compared with the 


traditional Back Propagation (BP) neural network, showing the feasibility and validity of the proposed approach. 
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1 Introduction 


In recent years, with decommissioning of a large 
number of nuclear facilities and increasing demand of 
waste treatments, the existed radioactive pipes need to 
be treated well. In order to sort and treat the alpha 
radioactive pipes before disposal, radioactivity survey 
is to be done to assess whether pipes meet limit 
requirements. Alpha particles released by the nuclides 
in the pipelines are absorbed partially by the pipe 
walls, which makes it difficult to detect the alpha 
pollution situation in the pipelines by external direct 
Alpha 
Detector (LRAD) technology provides an approach to 


testing methods. Recently, Long-Range 


get the alpha pollution information inside the pipelines. 


The LRAD detection system normally consists of five 
units, that is, sample detection (ion chamber and 
measurement chamber), air-driven part, power supply, 
signal acquisition and signal processing unit. This 
technique can detect indirectly alpha radioactivity by 


collecting ions inside pipelines that are produced by 
alpha particles. Thus, it can overcome defects of direct 
alpha detection such as short range and failing to 
penetrate facility walls. Many scholars have made 
preliminary studies to confirm that distance, length, 
diameter, radioactivity, wind speed and air flux 


1-5] Our statistical 


influence final measuring results! 
analysis to LRAD experimental results have shown 
that there is a nonlinearly relationship between testing 
parameters and measuring results. 

This paper puts forward a kind of improved 
with the uncertainty 


correction. The experiment data are collected as the 


genetic algorithm to deal 


following experiment conditions: the bias voltage of 
the current ionization chamber bias is 200 V. The 
diameters of carbon steel pipes are 43 mm, 48 mm, 
58mm, 66 mm, and 78 mm, and their are 10—160 cm 
which can be adjusted by rotary joints. The alpha 
sources are 24.05 Bq and 3200.00 Bq to establish the 
data training model and example detection model. The 
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results show that the improved genetic algorithm can 
effectively improve the precision of prediction. 


2 Methodology 


The GAANN model has three layers with m nodes in 
the input layer, A nodes in the hidden layer, and n 
nodes in the output layer. Firstly, the model is 
implemented in order to determine a basic state space 
of connection weights matrix. Secondly, the number of 
hidden nodes and connection weights matrix are 
encoded into a mixed string which consists of integer 
value and real value!”*!. In this paper the experimental 
data is divided into three parts: training sample 911, 912, 
cross-validation sample 921, g22 and testing sample 931, 
(32. Here we introduce our scheme: 

Step 1: Initialize connection weights which are 
within [—1,1] for training sample 911, gi2. Adjust the 
weights until the desired tolerance of error €11, €12 iS 
obtained. The maximum and minimum of weights are 


denoted as Umax and Umin, respectively. The value of 


weights are taken within [umin—O1, Umaxtd2], where 6), 


62 are adjustment parameters. 
minE Ay (t)-5,(0FP (1) 
i53 Yk POI <E; 


where i = 1, 2, 3 which corresponds to three data sets. 
P(t), y(t) are the desired output and real data. 

Step 2: Encode connection weights and number 
of hidden nodes. The hidden nodes are encoded as 
binary code string, 1 with connection to input and 
output nodes and 0 with no connection. The weights 
are encoded as float string, with string length H=mxh+ 
h+hxnt+n (m is the number of input nodes, n is the 
number of output nodes, h is the number of hidden 
nodes). Each string corresponds to a chromosome 
consisted of some gene sections, tabulated as follows: 

Step 3: Initialize a population of chromosomes. 
The length L of each chromosome equals to G+H, 
where G is the length of binary code of the number of 
hidden nodes and H is the length of real-valued code 
of connection weights. 


Table 1 Schematic diagram of encoding chromosome. Part A is encoded in binary type, and other parts in real value. These values 


change during training period 


l, l 0.2... 0.7 0.3, ...,0.1 0.2, ...,0.3 0.9,...,0.8 
A B C D E 
Step 4: Calculate fitness individually according ji ; 
X, =A eG (5) 


to Equation 2 below. 
F =1/(1+min £) (2) 
Step 5: Copy the highest fitness individual 
directly to a new offspring and select other individuals 
by the method of spinning the roulette wheel”! 

Step 6: Use basic crossover and mutation 
operations to the control code, namely, if a hidden 
node is deleted (added) according to mutation 
operation, the corresponding control code is encoded 0. 
The crossover and mutation operators of weights are 
encoded as follows: 

Crossover operation with probability pe 


Xi =c- X; +(l-¢,)-X; (3) 


i+] 


Xa = 1-6) tA (4) 


itl 


where X/, Xj,‘ are a pair of individuals before 
t+1 +1 i eal a 

crossover, X; , Xi" are a pair of individuals after 

crossover, c; is taken as random value within [0,1]. 


Mutation operation with probability pm 


where X/ is individual before mutation, X;"' is 
individual after mutation, c; is taken as random value 
within (umin 01 —X/ , Umax + 62 + X; ). 

Step 7: Generate the new population and 
replace the current population. The above procedures 
(step 4-7) are repeated until convergence conditions 
(min E:<epņ and min £3<¢;3) are satisfied, where A=1, 2, 
3 which corresponds to three data sets. 

Step 8: Decode the highest fitness individual, 
obtain corresponding number of the hidden nodes and 
connection weights, and output the prediction results. 


3 Model implementation and results 


In order to evaluate the performance of our model, we 
choose the representative 23Pu strong source (32.00 
Bq) and the weak source (24.05 Bq), respectively. And 
then establish a prediction model, in which the 
measured distance, tube length, diameter, wind speed 
and air flow construct the input layer, the ionization 
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voltage works as output layer. In this case, we choose 


ionization voltage value which can represent 


radioactivity intensity as a prediction of the output 
we find the linear 


value. In the experiments 


Table 2. Tabulated required model parameters 


relationship between radioactivity and ionization 
voltage value. After measuring the ionization voltage 
value, we take it for comparison with the standard 
source. And the radioactive intensity can be calculated. 


Name of variables Value Name of variables Value 
Training sample Q11 128 Error €;, of sample 911 0.05 
Validation sample 91 27 Error £12 of sample 917 0.10 
Testing sample 913 27 Error €)3 of sample 913 0.10 
Training sample 21 105 Error £2; of sample 921 0.15 
Validation sample 927 22 Error £2 of sample 27 0.30 
Testing sample 923 22 Error £23 of sample 93 0.30 
Number of input nodes 6 Crossover probability p, 0.50 
Number of hidden nodes 6 Mutation probability pm 0.10 
Number of output nodes 1 Population 20 
Learning rate 0.05 Iteration number 50 
Training epochs 1000 Training goal 10° 
200 
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© Training date 180 Validation R=0.9453 oe 
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Fig. 1 Measurement and prediction values under 3200.00 Bq. 


In the simulations a three-layered BP neural 
network is first employed to estimate basic state space 
of connection weights. The minimum and maximum 
values of weights are obtained, —1.21 and 0.96, 
respectively. Let 6;=—-0.09 and ô = 0.04, therefore, the 
range of weights is assumed to be within [—1.3, 1.0]. 
The number of input neurons is 5 and the hidden nodes 


200 60 80 100 120 140 160 180 200 


Measurement values 


are 6. The activation function adopted here from input 
to hidden layer is Sigmoid, while from hidden to 
output layer is Purelin function. For the proposed 
hybrid neural network, the following system 
parameters in Table 2 are applied to training samples 


and prediction. 


060201-3 


WANG Lei et al. / Nuclear Science and Techniques 24 (2013) 060201 


In order to test the model performance, we set 
training error ¢),=0.05, validation error &€ņ2=0.10, 
testing error €;3=0.10, under 3200.00 Bq alpha source. 
The training results are shown in Fig.1. 

As under the weak source environment, the 
measurement error is big. We set training error 
€91=0.15, validation error &€=0.30, testing error 
€23=0.30 under 24.05 Bq, as shown in Fig.2. 

For the purpose of comparison with other 
neural network models, such as the basic BP neural 


network, three types of errors, which are commonly 
found in many papers discussing these models, are 
also used here. Three types of errors are the mean 
absolute percentage error (MAPE), the maximum 
absolute percentage error (MAXAPE), the minimum 
absolute percentage error (MINAPE). 

As the initial connection weights and threshold 
have certain of randomness, so we have made two 
predicted models 50 times. 


2.4 
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Fig. 2 Measurement and prediction values under 24.05 Bq. 
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Table 3 Error comparisons of GA-based neural network and traditional neural 


Traditional BP network 


GAANN 
Confidence interval 0.2247-0.2816 
Numbers fallen in confidence interval 15 
Mean error 0.2373 
MAXAPE 0.6590 
MINAPE 0.0025 


0.2247-0.2527 
4 


0.2748 
0.7332 
0.0048 


Table 4 Error comparisons of GA-based neural network and traditional neural network in 3200.00 Bq environment 


GAANN 


Traditional BP network 


Confidence interval 


Number fallen in confidence interval 


Mean error 
MAXAPE 
MINAPE 


0.0667—0.0747 
12 

0.0707 

0.1555 

0.0098 
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0.0762~0.0845 


11 

0.0740 
0.1598 
0.0006 


WANG Liang et al. / Nuclear Science and Techniques 24 (2013) 060201 


We took its 95% confidence interval as neural 
network ensemble. The parameters list in Table 2. 
Characteristic parameters are tabulated in Table 3 and 
4. We can observe that the capability of approximation 
of this model is better than the traditional one for weak 
and strong radioactive sources. But in weak 
radioactive source, error precision is large. This is a 
shortcoming of the proposed model, which pushes us 
to present a better model and will report the results in 


future publication. 
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Fig.3 Measurement and prediction values under 24.05 Bq. 
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Fig.4 Measurement and prediction values under 3200.00 Bq. 


4 Conclusion 


In this paper, a GA-based neural network approach has 


been proposed for LRAD radioactivity survey research. 


Network structure is 
weights are adjusted through the implementation of 


optimized and connection 


genetic operators. The experiments with LRAD data 
have shown that the predictive performance of the 
proposed model is better than that of the traditional BP 
neural network. However, weak radioactive source 
exhibits the poor effect. The consideration of further 
improvements in model performance should include 
the following factors such as different time windows, 
different prediction horizon, crossover and mutation 
operators, classification of data sets, etc. This work is 
now under progress. 
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